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ABSTRACT 

The laminar incompressible viscous flow in a short cylindrical vortex chamber is 
analyzed numerically. The flow is assumed to enter the chamber priphery with an arbi- 
trary ratio of tangential to radial velocity components and to exit at a line sink at the 
chamber axis. The numerical technique permits consideration of both the tangential and 
radial momentum equations. Comparisons are made with solutions available for the 
limiting cases of zero swirl and infinite swirl. Results are presented for a wide variety 
of swirl ratios, Reynolds numbers, and chamber aspect ratios. The results include 
radial and tangential velocity profiles, pressure distributions, and points of radial flow 
reversal. 
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VISCOUS FLOW IN A SHORT CYLINDRICAL VORTEX 
CHAMBER WITH A FINITE SWIRL RATIO 
by Robert W. Horn beck 
Lewis Research Center 

SUMMARY 

The problem of laminar incompressible viscous flow in a short cylindrical vortex 
chamber is analyzed using a numerical technique. This technique permits the consid- 
eration of both the tangential and radial momentum equations as well as continuity, 
hence, the amount of swirl introduced into the flow at the periphery of the chamber may 
be varied. Velocity profiles and pressure distributions are obtained, assuming that the 
radial- and tangential -velocity distributions at the peripheral inlet to the chamber are 
uniform across the height. Values of Reynolds number Re based on the radial inlet 
velocity and chamber height between 20 and 2000 are considered. The swirl ratio S, 
defined as the ratio of tangential to radial inlet velocity, is varied from 0 to 50. Two 
values of chamber radius to chamber half -height (R^ = 5 and 10) are considered. 

Overshoots in the radial velocity within the end-wall boundary layers are found to 
occur somewhere in the chamber for all S > 1, and this overshoot (an increase in the 
radial velocity in the boundary layer over the velocity in the region far from the wall) 
is increased for increasing swirl and decreasing distance from the central axis of the 
chamber. For large values of Re. the end-wall boundary layers are very thin; hence 
the radial velocity distributions show sharp spikes in the boundary layer and are rela- 
tively uniform over the remaining channel height. For small Re, the end-wall boundary 
layers are sufficiently thick to fill the chamber height, except for the region close to the 
peripheral inlet. 

For s uff iciently high swirl ratios (greater than about 4 for most cases considered) a 
local reversal in the radial velocity will occur at some value of the radius, on the 
chamber centerplane. When this radial outflow occurs, the numerical scheme becomes 
unstable and must be terminated at the radius at which it occurs. 

The pressure distributions are found to be virtually identical to the potential flow 
distributions for S > 1, but are dependent on Re and R^ for S < 1. Even small 
amounts of swirl are found to have a relatively significant effect on the pressure drop 
between any two radial stations in the chamber. 



INTRODUCTION 


Considerable work has been done in recent years toward the description of flow in a 
confined vortex. A typical confined vortex configuration might consist of a cylindrical 
chamber with a fluid injected both radially and tangentially at the chamber periphery and 
withdrawn at the central axis. The work has been motivated at least partly by interest 
in a number of devices that utilize such a vortex configuration, including the Ranque- 
Hilsch tube, the magnetohydrodynamic vortex generator, the vortex nuclear reactor and, 
most recently, the vortex amplifier for fluidics applications. An excellent historical 
review of the vortex literature is given in Anderson (ref. 1). Reference 1 represents 
the most accurate solution to date for the velocity and temperature distributions inside a 
vortex chamber, but it is restricted to very high swirl so that the radial inflow velocity 
is negligible compared with the tangential (swirl) velocity in the stream between the end- 
wall boundary layers. Numerical solutions are then obtained assuming that the free 
stream is either irrotational or in solid -body rotation. 

The influence of various swirl ratios was investigated experimentally by Donaldson 
and Snedeker (ref. 2) but only for a cylindrical vortex chamber with a single end wall, 
rather than the fully confined vortex. In a recent analysis by Ostrach and Loper (ref. 3), 
finite swirl ratios were considered for the confined vortex, but it was necessary to 
linearize the basic equations in such a way as to ignore the effect of swirl ratio on the 
boundary -layer thickness. The present analysis retains the essential nonlinearities of 
the basic equations and accurately accounts for all of the effects of swirl ratio. 

Savino and Keshock have presented the most reliable experimental work to date on 
the radial and tangential velocity distributions in a confined vortex (ref. 4). Although 
the flow was turbulent, this investigation allows at least a qualitative comparison with 
the work of Anderson (ref. 1) and the present study for the case of high swirl ratios 
015). 

In the present investigation, viscous incompressible flow in a confined vortex is ana- 
lyzed by using a numerical marching technique that permits full consideration of both tan- 
gential and radial momentum at all points in the chamber, thus allowing consideration of 
any swirl ratio, as well as viscous effects over the entire chamber. In order to employ a 
numerical marching technique, it is necessary to drop certain terms in the basic equations 
This restricts the solution to short vortex chambers (chambers in which the radius is larg- 
er than the height) . One drawback to the use of this numerical scheme is that the numer- 
ical technique becomes unstable when radial outflow is present. As a result of this, the 
solution may only be carried inward from the periphery of the chamber until radial out- 
flow begins to occur. Despite this shortcoming, the method does permit the detailed study 
of a large number of interesting cases, particularly since the real area of interest in this 
investigation is for finite swirl ratios, where radial outflow is less likely to occur. 
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PROBLEM AND BASIC EQUATIONS 


The problem configuration for a typical short vortex chamber is shown in figure 1. 
The fluid is assumed to enter the chamber both radially and tangentially with velocity 
distributions that are uniform over the height of the chamber. The flow is assumed to 
be laminar and incompressible, with constant properties. The fluid is removed at the 
center of the channel by what is assumed to be a line sink. In the actual chamber, 
central exit tubes of relatively small diameter must be provided at the top and bottom 
of the chamber. These, presumably, will not have a large effect on the flow pattern in 
the chamber except close to the center. Axial symmetry is assumed throughout. 

The basic equations of motion for the steady flow of an incompressible, constant 
property fluid in cylindrical coordinates are 
r -direction momentum: 


p L *i + u 

\ 3z 3r 



dr \3z 2 9r 2 r9r r 2 / 


( 1 ) 


0-direction momentum: 


P 




3z 


3r 


3 2 v + 3 2 v + 1 3v _ v \ 
,3z 2 3r 2 r9r r 2 / 


( 2 ) 


continuity: 


3w h 1 3(ur) = 0 
3z r 3r 


(3) 


9 o 

When an order of magnitude analysis is made and terms of € (h /r|) are neglected, 
the basic equations are reduced to 
r-direction momentum: 0 


P 


iw — + u — 


0Z 


3r 




(4) 
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9 -direction momentum: 


P 




continuity: 


3w + 1 3(ur) _ 0 
3z r 3r 

Equations (4) to (6) apply only for small (h/r^, that is, short vortex chambers. The 
boundary conditions on equations (4) to (6) are 

u(r z) = -Uj = Constant 

u(r, h) = 0 

— (r, 0) = 0 
3z 

v(r p z) = Vj = Constant L 

v(r,h) = 0 
— (r, 0) = 0 

3z 

w(r,h) = 0 
w(r; 0) = 0 

In order to attempt a numerical solution, it is first necessary to put the basic 
equations into a dimensionless form which reduces the number of parameters that 
appear to a minimum. The dimensionless variables chosen are 
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U 1 


V = — 

U 1 


W = — 

U 1 


p_ p - p l 

P uf 


( 8 ) 


Z 


z 

h 



Inserting the dimensionless variables (8) into equations (4) to (6) results in the following 
dimensionless forms of the basic equations: 


W 3U , U 3U V 2 - dP , 1 d 2 U 

3Z 3R R dR Re az 2 

w av + u 9V UV - 1 9 2 V 

az 3R R Re az 2 

aw ^ i a(UR) = 0 

az R 3R 


( 9 ) 

( 10 ) 

( 11 ) 


where Re = pujh/ p is a Reynolds number based on the channel height. This parameter 
is a natural result of the dimensionless variables chosen earlier. The boundary con- 
ditions (7) in dimensionless form become 
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U(R X ,Z) = -1 


*•% 


U(R, 1) = 0 
— (R, 0) = 0 

az 

V(R 1 ,Z)=Ii = S 
U 1 


where 


V(R, 1) = 0 f 


— (R, 0) = 0 

az 


( 12 ) 


W(R, 1) = 0 
W(R, 0) = 0 


R 


1 


^1 

h 


The parameters (including boundary conditions) that must be considered for this problem 
are a geometry parameter r^/h, the swirl ratio v^/u^, and the Reynolds number 
pujh/ju. 


Finite Difference Representation 

The finite difference grid employed in this analysis is shown in figure 2. Note that 
increasing j subscript corresponds to decreasing radius, since the marching direction 
will be from the outside of the chamber toward the center. Equations (9) to (11) will now 
be expressed in an implicit finite difference representation at a point (j + l,k): 
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w. U j+l,k+l U j+l,k-l u U j,k " U j+l,k _ V j,k V j+l,k 

3,k 2 AZ 3 ’ k AR R j + 1 


P J - P i + 1 . 1 

AR Re 


U j+l,k+l " 2U j+l,k +U j+l,k-l 


(AZ)- 


(13) 


W. 


V, 


i»k 


j+l, k+1 ~ V j+l,k-l | n V j,k ~ V j+l,k + V j+l,k U j,k 


2 AZ 


1,k 


AR 


R j+1 


_1_ 

Re 


V j+l,k+l ~ 2V j+l,k + V j+l,k-l 


(AZ)- 


(14) 


W j+l,k+l W j+l,k + 
AZ 


R j+1 


U. , R. - U- , , R . . 

J, k j J+l, k j+l 


AR 


= 0 


(15) 


Equations (13) to (15), written for k = 0(l)n (i.e., for k ranging from 0 to n in steps 
of 1), constitute a complete set of (3n + 3) equations in the (3n + 3) unknowns 

U j+1,0’ • * ’ U j+1, n’ V j+1,0’ * ‘ * V j+1, n’ W j+l,l’ * ' • W j+l,n ; and P j+r 


Method of Solution 


In order to obtain a solution to equations (13) to (15), the first step will be to re- 
arrange equations (13) and (14) into a more useful form. Equation (13) becomes 


K k 

1 

u. , , . + 

r-°i,k , 2 1 

TJ. „ , + 

[ W 3, k 1 1 

2(AZ) 

Re (A Z ) 2 _ 

u ]+l,k-l + 

_ AR Re(AZ) 2 _ 

3+1, k + 

_ 2 ( AZ) Re(AZ) 2 


U, 


j + l, k+1 
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Equation (14) becomes 


Fv 

1 

V j+l,k-l + 

F'°i. k , D J. k + ^ ] 

2(AZ) 

Re(AZ) 2 _ 

_ AR R j+1 Re(AZ) 2 _ 


'j+l,k 


W ; 


j, k _ 


V« 


2(AZ) 


Re(AZ)^ 


j+l,k+l 


- u i.k v i, k 

AR 


(17) 


Equations 


U j+1,0’ ' 
continuity 

additional 


(16) and (17) for k = 0(l)n comprise (2n + 2) equations in the (2n + 3) unknowns 
' • U j+1, n’ V j+1,0’ ' ■ * V j+1, n’ and P j+1- Rather than employing the 
difference equation (15) to complete the set thus introducing the W^ +1 k ’s as. 
unknowns, it is advantageous to introduce the integral form of continuity 


2?rR j+l 




27tRj 


Equation (18) may be written in finite difference form as 




(19) 


Equation (19) together with equations (16) and (17) written for k = 0(l)n now completes 

the set of (2n + 3) equations in the (2n + 3) unknowns U. +1 ^ Q , . . . Uj + i n >* V j+1, O’ ' ‘ ‘ 

V- , : and P. , . This is a significant reduction in number from the (3n + 3) equations 

3 + 1 , n 7 ] +1 

which it would be necessary to solve if the continuity difference equation (15) were used. 
The use of an integral continuity equation of this type in finite difference analyses of 
confined flow problems was apparently first used by Bodoia and Osterle (ref. 5), and 
the techniques used herein have their foundations in that work and in the work of 
Hornbeck, Rouleau, and Osterle (ref, 6) in which two momentum equations and continuity 
were solved simultaneously for the entrance flow in a porous tube. 

The complete set of equations may be written in matrix form as 
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^0 

y 0 


^0 
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U j+1, 0 


Vo 

“1 

h 

y l 


*1 


V 


U J + 1, 1 


n 


“2 

^2 y 2 


$2 


V 


U j+1, 2 


V 2 



. 


- 


- 


- 


- 



- - - 



- 

- 


- 


- 



“n-1 ^n-1 y n-1 



*n-l 

V 


U j+l,n-l 


Vn-l 



a n 




V 


U j+1, n 


V n 




B 0 

G 0 


0 

X 

v j+i, o 

= 

F 0 




A 1 

B 1 G 1 


0 


V j+1, 1 


F 1 





A 2 B 2 

G 2 

0 


V j+1, 2 


F 2 






_ 

_ 


- 


- 






A n-1 B n-l G n-l 

0 


V j+l,n-l 


F n-1 






A n B n 

0 


V j+1, n 


F n 

1 

2 

2 - - 2 2 

0 

0 0 0 

0 0 0 

0 


P j+1 

1 

: 

Q 


where 




AR 


Re(AZ)^ 


^0 


-2 


Re(AZ)^ 


V s - 


-1 

AR 


"k = 


- w i.k 

2(AZ) 


Re(AZ) / 




“ U j,k 

= — Al — + 

AR 


Re(AZ^ 


M 


(k*0) 
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I'k 


= . 
2(AZ) 


(k#0) 


Re(AZ)^ 


h = - 


V j,k 

R i+ i 


r k = 


-U 2 , P. 

3,k- ] 

AR 


b .Ao/U 2 

0 


AR R j+1 Re(AZ) 2 


G 0 = 


Re(AZ)^ 


W i k 

A _ _ 

k " 2(AZ) 


Re(AZr 


(k^O) 


B k = 


V + V + _j^ 

AR R j+1 Re(AZ) 2 


(k#0) 


W i k 

G k = — ~ — “ 
K 2(AZ) 


Re(AZr 


(k^O) 


F = - U 3> kV J> k 
k AR 


2R-, 


Q = 


(A Z )R • 


3+1 


The symmetry conditions have been included in equation (20) for k = 0. 
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The solution to equation (20) may be obtained by any of the standard methods used 
for simultaneous linear equation solutions, such as Gaussian elimination, a modification 
of Gauss -Seidel iteration, or inverting the coefficient matrix and multiplying the right 
hand side column vector by the inverse. 

Once the solution to the set of equations has been obtained, it will be necessary to 
obtain the values W. +1 j, . . . W i+1 n before proceeding. This may be readily 
accomplished by solving equation (15) for Wj +1 k+ p yielding 

w i + i,k + i • w i+ i,k - £5^ [ u j,k R r u j+i,k R j+i] (21) 

This solution for w j + j ^+1 ma ^ be sim P 1 y advanced in the direction of increasing k, 
starting from k = 0 and proceeding step by step. 

After all values of Uj + ^ Vj + j ^j+1 k’ ^j+l ^ ave ^> een found, another 
step AR inward along the radius may be taken and the entire process repeated. 

Because of the computational difficulties that occur when the boundary layers are 
very thin and the velocity gradients high, it is advantageous to use a considerably smaller 
AZ mesh close to the wall than that used toward the midplane of the chamber. The 
details of the method used for changing the mesh size are given in references 7 and 8. 

RESULTS AND DISCUSSION 

The equations derived in the previous section were solved on a high-speed digital 
computer. The axial mesh sizes used were AZ = 0. 1 toward the chamber centerplane 
and AZ = 0. 015 in the region close to the walls, that is, Z > 0. 7. For a few of the 
high Reynolds number cases where the boundary layers were very thin, it was necessary 
to, instead, use AZ = 0. 005 for Z > 0. 9 and AZ = 0. 1 for Z < 0. 9. The radial 
mesh sizes used were AR = 0. 025 for the R^ = 5 cases and AR = 0. 05 for the 
R^ = 10 cases. 

The results presented include radial- and tangential -velocity profiles and pressure 
distributions. The values of parameters considered include R^ = 5 and 10; Re = 15, 

20, 50, 200, and 2000; and S = 0. 2, 0.5, 0.75, 1, 2, 4, 10, and 50. A rather detailed 
presentation of radial-velocity distributions is given for S ^ 1, but only a few repre- 
sentative tangential-velocity profiles are shown because of the general similarity of 
these profiles. Pressure distributions are illustrated for the full range of parameters. 
Also included is a graph of the radii at which radial flow reversal takes place as a 
function of swirl ratio and Reynolds number, and a set of comparisons with existing 
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solutions. 

Figure 3 compares the profiles of the radial velocity component obtained in this 
investigation with those obtained by Anderson (ref. 1) for very high swirl. Anderson’s 
profiles were derived for two types of imposed rotational motion in the flow region 
outside the end-wall boundary layer: (1) irrotational vortex motion and (2) solid-body 
rotation. Anderson’s profiles revealed a tendency for a small magnitude flow reversal 
to develop only when the main stream was in a solid-body rotational mode. This re- 
versal occurred near the outside diameter of the chamber in the region 0. 8 < R/Rj < 

1. 0. With the irrotational vortex motion case, Anderson found no tendency for the 
development of flow reversal. 

Typical radial velocity profiles obtained by Anderson and from the present numerical 
solution are compared in highly magnified form in figure 3(a) for Re = 15 and S =20. 
For these conditions, the irrotational and solid-body rotation solutions of Anderson are 
quite similar, and the solution obtained in the present investigation lies between the two 
and favorably agrees with both of the profiles in general shape and location of the peak. 
The agreement is best with the irrotational vortex-solution profile. The low Reynolds 
number solutions obtained in the present investigation behave in a qualitative way like 
the solid-body rotation solutions of Anderson, in that they show a distinct tendency to 
flow reversal quite close to the chamber periphery. The value of S = 20, which was 
used to obtain the results shown in figure 3(a) for the present investigation, was the 
largest value of swirl ratio that could be used for Re = 15 and still allow the radial 
position R = 9 to be reached without reversal. 

Figure 3(b) shows a comparison of the solution obtained from the present investi- 
gation with the solutions of Anderson for Re = 2000 and S = 50. The profiles are 
very highly magnified and the agreement of the present solution with the irrotational 
vortex solution is excellent. The small discrepancies that are present are caused at 
least partly by the fact that there is always some radial forced flow in the present 
investigation, whereas any radial flow present in Anderson's solution is entirely induced 
flow. This results in a somewhat higher radial flow rate for the results of the present 
investigation. This discrepancy would tend to decrease with increasing swirl ratio. 

The solid-body rotation solution does not give good agreement with the present solution. 
The profile shown for the solid-body rotation solution has reversed at the chamber 
centerplane, but the profile from the present solution shows no tendency toward flow 
reversal at this high Reynolds number and differs significantly in the shape and magnitude 
of the peak from the solid-body rotation solution. 

The present investigation thus yields solutions that fall somewhere between the 
irrotational vortex and solid -body rotation solutions of Anderson, with some character- 
istics of both. The best agreement over the entire range of parameters appears to be 
with the irrotational vortex solution. Further evidence that the vortex flow outside of 
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the end-wall boundary layers is essentially irrotational may be found later in this 
report in the discussion of tangential -velocity profiles. 

Figure 4 compares results of the present investigation and those of two other inves- 
tigations for the radial pressure distribution with no swirl, that is, radial convergent 
flow between parallel plates 0 The solutions used for comparison are those of Moller 
(ref. 9) and Comolet (ref. 10). Results identical with those of Comolet were also found 
independtly by Levesey (ref. 11). These analyses are momentum integral solutions in 
which some form of parabolic velocity profile is assumed. The agreement of the present 
numerical solution with the integral solutions is quite good. The numerical solution 
agrees very well with that of Moller for R > 3 and with that of Comolet for R < 3, 

In the region R > 5, the pressure drop of the present solution is higher than that 
obtained from either of the integral solutions, which is due at least in part to the higher 
pressure drop that results from the uniform inlet velocity profile used in the present 
analysis as opposed to the parabolic profile employed in the other investigations. 

A solution including inertia for zero swirl flow between parallel plates was carried 
out by Peube (ref. 12) using a series expansion, but, for the values of parameters 
considered herein, the truncated form of the expansion given by Peube is not adequate, 
and presumably many more terms would be required for reasonable accuracy. No 
comparison, therefore, will be made with this solution. 

Figure 5 shows radial velocity distributions for Reynolds numbers of 50, 200, and 
2000. For each Reynolds number, values of R^ of 5 and 10 are considered. The swirl 
ratio S appears on each graph as a parameter that varies from 1 to 50 for Re = 50, 
from 1 to 10 for Re = 200, and from 1 to 4 for Re = 2000. Larger values for the upper 
limit of S for Re = 200 and Re = 2000 were not felt to be practical since the corre- 
sponding Reynolds number based on the tangential velocity would exceed any reasonable 
limit for laminar flow. 

Figure 5(a) illustrates the radial -velocity overshoot in the boundary layer, which is 
the most striking characteristic of vortex flow. The higher the swirl ratio, the greater 
the tendency toward this overshoot, while the radial velocity near the chamber center- 
plane is correspondingly decreased. This effect is also accentuated as the axis of the 
chamber is approached as may be seen by comparing figures 5(a-l) and (a-2). This 
overshoot phenomenon is caused by the tangential velocity decreasing rapidly in the 
boundary -layer region, which reduces the centrifugal effect in that region and allows the 
radial pressure gradient to drive the flow inward along the chamber wall. 

The end-wall boundary -layer thickness is an inverse function of the Reynolds 
number, and, for Re = 50 and R-^ = 10 in (fig. 5(a-l), the boundary layers have 
reached almost halfway to the chamber centerplane at only one -tenth of the distance from 
the periphery to the axis of the chamber. 
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One significant result of the swirl is that, for even moderately high swirl ratios 
(S > 4), a significant portion of the radial flow is carried in the boundary layers close 
to the end walls. For example, in figure 5(a-2) we find that for S = 4 about one-third 
of the total radial flow is confined to the region between Z = 0. 8 and the end wall when 
a radial position R = 5 halfway to the chamber axis has been reached. 

Figures 5(a-3) and (a-4) show results for the same Reynolds number of 50 used to 
obtain the results shown in figures 5(a-l) and (a-2) but for a chamber with a smaller 
dimensionless radius Rj = 5. The boundary layers occupy a slightly smaller fraction 
of the total chamber height for the R^ = 5 cases than for the R^ = 10 cases at the 
same relative radial location R/R^ when all other parameters are held constant. This 
seems to be the only significant effect of varying R 1 . 

Figure 5(b) shows a set of radial -velocity profiles with Re = 200 and Rj = 10 and 
5. These results differ from those for the Re = 50 cases, primarily in the effect of 
boundary -layer thickness on the velocity profile. Since the boundary layers become 
thinner for higher Reynolds numbers, the velocity overshoot is confined to a smaller 
region close to the end walls. This thinner boundary layer has less effect on the 
centerplane radial velocity and thus there is less tendency toward flow reversal at the 
higher Reynolds numbers. The thinner boundary layer for Re = 200 causes a con- 
siderable difference in appearance between the Re = 50 and the Re = 200 results, 
but for S = 4 (fig. 5(b-2)), about one-third of the radial flow is confined to the region 
between Z = 0. 8 and the end wall, which was the same result found for the correspond- 
ing Re = 50 case. As before, the only effect of varying R^ seems to be to cause a 
slightly thinner boundary layer for the R^ = 5 cases shown in figures 5(b-3) and (b-4) 
as compared with the R 1 = 10 cases shown in figures 5(b-l) and (b-2). 

The radial velocity distributions for Re = 2000 shown in figure 5(c) continue the 
same trends discussed previously for the lower Reynolds numbers. For the case shown 
in figure 5(c-2) again about one-third of the radial flow rate is contained between 
Z = 0. 8 and the chamber end wall for R = 5 and S = 4. 

Some representative tangential -velocity profiles are shown in figure 6. Figure 6(a) 
illustrates a low Reynolds number case which results in thick boundary layers. Despite 
this thick boundary layer, the centerplane tangential velocity follows a l/R variation 
quite closely over the entire range of R. This indicates that the flow is close to 
irrotational at the chamber centerplane. For the Re = 2000 case shown in figure 6(b), 
the boundary layer is much thinner, and the flow outside the boundary layer is essenti- 
ally irrotational. Tangential -velocity profiles for the other values of the parameters 
have this same general behavior, with no velocity overshoots or inflections present. 

Figure 7 illustrates the relative radial position R/R 1 at which radial -velocity 
reversals at the chamber centerplane begin for Rj = 10. The point at which this 
reversal occurs is a function of Re and S. With increasing swirl ratio, the radial - 
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velocity reversal takes place closer to the chamber periphery. The reversal is retarded 
by higher Reynolds numbers due to the resultant thinner boundary layers. Reversal 
begins to take place somewhere in the chamber at S = 4 for Re = 50 and at S = 20 for 
Re = 200, while for Re = 2000 no reversal occurs, even for S = 50. Results for R^ = 5 
would appear much the same as those shown in figure 7, but all points would be lowered 
slightly since flow reversal would have less tendency to occur due to the slightly thinner 
boundary layers that result for the R^ = 5 cases. 

The pressure distribution results are shown in figures 8 and 9. Figure 8 shows 
dimensionless pressure as a function of relative radial position R/R^ for swirl ratios 
greater than 1. For this range of swirl ratios, the predominant factor in the pressure 
drop is the centrifugal effect, with the frictional pressure drop providing negligible 
contribution. As a result of this, the single set of curves shown in figure 8 is sufficient 
to describe the pressure distribution completely independent of Re and R^. The 
pressure varies strongly along the radius for these swirl ratios as the tangential velocity 
increases at the expense of the static pressure. The pressure drop becomes essentially 
infinite as the center of the chamber is approached. 

Figure 9 shows the radial pressure distributions for swirl ratios less than 1. In 

this range, the pressure distribution is a composite of frictional and centrifugal effects. 
Results are presented for three Reynolds numbers: Re = 20, 200, and 2000. The 
results for Re = 20 shown in figures 9(a) and (b) are somewhat different in shape for 
Rj^ = 5 and R 1 = 10. For the Re = 200 and Re = 2000 cases shown in figures 9(c) 
and (d), the same set of curves results for both R^ = 5 and R^ = 10, indicating that 
R^ dependence is important only for small Re. Therefore, in figures 9(c) and (d), the 
pressure distributions have been plotted as a function of R/R^. To give some idea of 
the effect of adding swirl to the flow, the increase in S from 0 to 1 gives an increase 
in the pressure drop at R/R^ = 0. 2 ranging from about 25 percent for the Re = 20, 

R^ = 10 case shown in figure 9(b) to almost 100 percent for the Re = 2000 case shown 
in figure 9(d). 

Since the experimental work of Savino and Keshock (ref. 4) was for a turbulent flow 
situation, no direct comparisons with the present analysis can be made. Qualitative 
comparisons can be made however, and the experimental results for the radial velocity 
distributions show the same tendency toward a high radial velocity overshoot close to 
the end wall which increases in magnitude with decreasing radius. These experimental 
results were obtained for a swirl ratio of about 15. All the experimental profiles show 
small radial outflow at or near the chamber centerplane, and the results obtained in the 
present investigation also show a tendency to flow reversal, particularly for high swirl 
ratios. The tangential velocities near the chamber centerplane in both the experimental 
work and the present analysis increase with decreasing radius with close to a l/R 
variation. This variation is essentially that of an irrotational vortex. The depression 
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in the tangential -velocity profiles near the chamber centerplane, which was noted in the 
experimental work, does not occur in the present analysis since this effect occurs only 
in the presence of significant radial outflow. 


CONCLUSIONS 

The laminar, incompressible, viscous flow in a short cylindrical vortex chamber 
was analyzed numerically. The flow was assumed to enter the chamber with an arbi- 
trary swirl ratio and to exit through a line sink at the chamber axis. The results may 
be summarized as follows: 

1. A velocity overshoot occurred in the radial velocity distribution close to the end 
wall of the chamber. This overshoot increased with increasing swirl ratio and decreasing 
radius. 

2. The radial velocity at the chamber centerplane decreased with decreasing radius 
and tended to become negative. This effect was strengthened with increasing swirl 
ratio. 

3. The boundary -layer thickness on the chamber end walls increased with decreasing 
Reynolds number and tended to weakly increase with increasing dimensionless chamber 
radius. 

4. The tangential flow pattern outside the end-wall boundary layers was essentially 
an ir rotational vortex. 

5. The radial pressure variation could be adequately found by assuming the flow 

to be frictionless and the vortex to be irrotational as long as the swirl ratio was greater 
than 1. 

6. For swirl ratios less than or equal to 1, both swirl and viscous effects made 
significant contributions to the radial pressure drop. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 27, 1968, 

129-01-05-20-22. 
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APPENDIX -SYMBOLS 


A,B,F,G 

h 

P 

P 

Pi 

Q 

R 

Re 


S 

U 

u 

U 1 

V 

v 

V I 
W 

w 

Z 

z 

P,Y 

AR 

AZ 

<P 

P 


coefficients defined in matrix eq. (20) 
chamber half -height 
dimensionless pressure, (p - p^)/pu^ 
pressure 

pressure at peripheral inlet of chamber 

constant defined in eq. (20) 

dimensionless radius, r/h 

Reynolds number, pu^h/p 

dimensionless outer radius of chamber, r^/h 

radial coordinate 

radius of chamber 

swirl ratio, v-j/u^ 

dimensionless radial velocity, u/u^ 

radial velocity component 

radial inlet velocity at chamber periphery 

dimensionless tangential velocity, v/u^ 

tangential velocity component 

tangential inlet velocity at chamber periphery 

dimensionless axial velocity, w/u^ 

axial velocity component 

dimensionless axial coordinate, z/h 

axial coordinate 

coefficients defined in eq. (20) 

radial finite difference mesh 

axial finite difference mesh 

coefficient defined in eq. (20) 

density 

viscosity 



Subscripts: 

j denotes radial (R) location in finite difference grid; those quantities having subscript 
j are known, those with subscript j + 1 are unknown 

k denotes axial (Z) location in finite difference grid; k = 0 is chamber centerplane 

n denotes axial location one step (AZ) removed from channel end wall; channel end 

wall is at k = n + 1 
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Axial coordinate, Z, dimensionless 




(a) Concluded. 


Figure 5. - Continued. 
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Figure 8. - Radial-pressure distributions for high swirl. 


Swirl ratio, S 

Figure 7. - Radial flow reversal points. Dimensionless chamber radius, 10. 


Relative radial location, R/R 






(c) Reynolds number, 200; dimensionless chamber radius, 10 and 5. 



Relative radial position, R/R x , dimensionless 
(d) Reynolds number, 2000; dimensionless chamber radius, 10 and 5. 
Figure 9. - Concluded. 
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